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ABSTRACT 

The Taiwanese-American Occultation Survey (TAOS) operates four fully au- 
tomatic telescopes to search for occultations of stars by Kuiper Belt Objects. It 
is a versatile facility that is also useful for the study of initial optical GRB after- 
glows. This paper provides a detailed description of the TAOS multi-telescope 
system, control software, and high-speed imaging. 
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1. Introduction 



The Taiwanese American Occultation Survey (TAOS) continuously monitors up to 
~1,000 stars with four telescopes at 5 Hz, for the purpose of detecting occultations of these 
stars by small (>1 km) objects in the Kuiper Belt and beyond. Installation of the first three 
telescopes and the first version of the control software were completed in January 2005, and 
the survey was officially started shortly thereafter. Installation of the fourth telescope was 
completed in January 2007, and four telescope operations commenced in August 2008. The 
TAOS Collaboration has finished analysis of the first two years of three-telescope data, which 
co mprises over 150,000 star-hours of observations. The results of this analysis are presented 



in (IZhang et al.ll2008a[ l. 



The principal science questions addressed with this survey are: 



• What are the number and size distribution of small (~0.5-10 km in diameter) bodies 
in the trans-Neptunian region? 

• Is there an extension of the Kuiper Belt beyond 50 AU comprising bodies too small to 
have been detected in direct surveys? 

• What are the number and size distribution of objects (e.g. Sedna) in the Extended 
Disk of the Solar System? 



The signatures we are searching for are brief reductions in flux from background stars due 
to occultations by small Kuiper Belt Objects (KBOs) and Extended Disk Objects (EDOs). 
The occultation survey technique is uniquely suited to the detection of such objects, since it is 
capable of det ecting objects that are much fainter than the magnitude limit of any plausible 
direct survey (iBaileyi Il976l ). However, the occultation technique is difficult to implement 
in practice, so the project includes some important technology development goals. For 
mo re details on the o cculta ti on techniqu e appli e d to sear c hes fo r outer Solar System object s 



see 



Bickerton et al. (2008). Nihei et al. 



2007h. ICooravl J2003h. ICoorav fc Fari^ jHoJ, 



Roques et al.l (l2006l ). IChang et al.l f|2007f). a nd iRoaues et al.l (20031 ). For an overview o f 



Alcock et al 



the TAOS Project, see iLehner et al.l (120061 ). IChen et all (120071). an d 
For more information on th e Kuiper Belt , see |j"ewitt &: Luij (119931 ) , ILuu fc JewittI _ _ . . , 
Bernstein et al.l (120041 ) , and iTrujillo et al. J200lh . For more information on the Extended 



( 2003 
J2002f ) 
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Disk, see iBrown et al.l (120041 ). iKenyon fc Bromleyl (12004 ). iMorbidelli fc LevisonI (120041 ). and 



Brasser et all (120071 ). 



In addition to the primary science goals stated above, a rapid response capability was 
implemented on the TAOS system to observe optical counterparts to Gamma Ray Burst 
(GRB) events reported to the ORB Coordinates Networlfl (GCN) The TAOS syste m ob- 
served its first optical counterpart to a GRB in October of 2007 (IWang et al.ll2008l ). The 
follow-up capability will be described in §5.7[ 



1.1. Design Criteria 

Implementation of an occultation survey for detection of small KBOs presents a num- 
ber of challenges. First, such events have very short durations, typically about 200 ms 
( iNihei et al.l 120071 ). Second, such events are extremely rare. Estimates of occultation 
event rates are model dependent and range from as low as 10~^ events per star p er year 
(IBernstein et al.ll2004l ) to as high as 10~^ events per star per year (jPan fc Sarill2005l ). (Note 



that when this survey was designed in 1996, occultation event rates as high as 10~^ per 
star per year were pr edicted for objects down to 1 km in diameter, but the results of 
Bernstein et al.l (120041 ) showed a likely break in the size spectrum at diameters of about 
30 km, which lowered the expected event rate.) 

The TAOS system was thus designed to meet the following criteria: 

• Make photometric measurements at a high sampling rate (5 Hz). 

• Continuously monitor enough stars (~1,000) to obtain a significant event rate. 

• Low (<0.1 per year) false positive rate. 

The rapid photometry is achieved by means of the innovative use of an otherwise conventional 
CCD camera (which will be described in §1]). An adequate number of stars can be imaged 
by a modern CCD camera mounted on a small (50 cm), fast (F/1.9) telescope with a wide 
field of view (3n°). 

Achieving the goal of a low false positive rate was a primary driver of the design. 
Our goal was to be able to measure a rate of stellar occultations by KBOs of 1 event per 
1,000 stars per year of observing time. We make ~10^° photometric measurements per year. 



^http : //gen . gsf c . nasa . gov/| see also iBarthelmv et al.1 (|1998[ ) 
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from which we want to derive less than one spurious occultation. The false alarm probability 
per observation must therefore be exceedingly small if the results of the survey are to be 
interpreted with confidence. We control false alarms by requiring joint detections in an array 
of four independent robotic telescopes. False positives of terrestrial or atmospheric origin 
(such as birds, aircraft, extreme atmospheric scintillation events, etc.) are removed by this 
requirement. False positives due to main belt asteroids can be eliminated by automated 
follow-up of detected events at a nearby, larger telescope. However, this would require 
a real-time analysis pipeline, which has not yet been implemented for this project. The 
telescopes are spaced far enough apart that scintillation events are unlikely to affect all of 
the telecopes simultaneously. Note, however, that the telescopes a re close eno u gh th at any 
occultation events would be detected by all four telescopes (see iNihei et al.l (120071 ) for a 
detailed description of occultation events from objects in the Outer Solar System). Finally, 
note that the redundancy in the multi-telescope system also greatly enhances the credibility 
of the reported results. It would be difficult to credibly report on rare occultation events 
that were seen in only a single telescope. 

This paper will specifically address the first two design criteria listed above. A detailed 
description of the survey system (both hardware and software) will be presented, and the 
implementation of the survey system to address the design criteria will be discussed. We have 
developed statistical analysis techniques which make efficient use of the multi-telescope data 
to limit the false positive rate and provide a credible result, however, discussion of these 
techniques is outside the scope of this paper. An introduction to th e statistical analysis 
methods used by the TAOS project can be found in iLiang et al.l (120041). and a discus sion of 
the application of these techniques will appear in a future paper (ILehner et al.ll2009l ). 



2. The Site 

The telescopes are installed atop Lu-Lin Mountain (longitude 120° 50' 28" E; latitude 
23° 30' N, elevation 2850 meters), in the Yu Shan (Jade Mountain) area of Taiwan. The 
median seeing is 1.3", with ~100 clear nights per year. Most of the poor nights occur during 
rainy season (May through September). Since the beginning of 2005, TAOS has averaged 
~400 hours per year of actual observing, an average of four hours per clear night. High 
humidity is typically the reason the observing hours are limited on clear nights. (Plans are 
currently being developed to start a new survey at a drier site.) 

The telescopes have separations ranging from 6 meters to 60 meters (see Figure[T]). These 
separations are sufficient to eliminate most, if not all, false positives arising from extreme 
scintillation events. Note that the original design for the TAOS project had the fourth 



telescope sited at a different peak located 7 km to the west of the Lu-Lin site ( lKing||200ll ). 
The reason for this was that differences in the timing of an event would give a measurement of 
the relative velocity of the occulting object, thus allowing the object distance to be estimated. 
This plan was abandoned for several reasons. First, the development of the second site was 
logistically very difficult, and hence extremely expensive. Second, an occultation shadow 
crossing the Lu-Lin site would cross the second site only during part of the year, due to the 
obliquity of the Earth. Third, our 5 Hz sampling rate is not nearly fast enough to give a 
velocity measurement sufficiently accurate to provide a reasonable estimate of the distance. 
Finally, it was determined that a four-telescope system would provide significantly better 
statistical identification of false positives than a thre e-telescope sy s tem, as well as enable 
detection of smaller objects. This will be described in lLehner et al.l (120091 ). 



3. The Hardware 

3.1. The Telescopes 

The four identical telescopes were manufactured by Torus Technologies (now Optical 
Mechanics, Inc.). The TAOS telescopes were designed to be compact, and the optical design 
was based on the system requirement to cover a large number of stars using a camera with 
a single CCD chip. Considering the affordable cost and ease of manufacturing, a 50 cm 
aperture was chosen. Each telescope is a Cassegrain system with the focal plane located 
behind the primary mirror. In order to provide the fast converging beam for a wide field of 
view, a parabolic primary mirror with F/1.5 is used. Combined with a large F/3.6 spherical 
secondary mirror (25 cm diameter), the overall F ratio is about 1.9 with an effective field 
of view of 3n°. With a 14 fim pixel CCD, the plate scale is about 3 arc sec/pixel, which 
undersamples the seeing of the site. The effective collecting area is about 1875 cm^.The wide 
field correctors consist of five lenses to provide a fiat focal plane. Four different materials 
were used for the chromatic correction. The design criterion was to provide 60% of the 
enclosed energy within one pixel area within the passband between 500 to 750 nm. The 
design and modeling of the telescope was done using the ZEMAX^^software packag(||. The 
optical layout is shown in FigureOand the detailed parameters for the design are summarized 
in Table m 



Table 1. Optical components of the TAOS telescopes 



Element 


Radius of Curvature (mm) 


Axial Thickness (mm) 


Clear Aperture (mm) 


Conic Constant 


Material 


Primary mirror 


-1496.332 


-418.142 


505.910 


-1 


Pyrex 
Air 


Secondary mirror 


-1795.977 


304.360 


244.953 





Pyrex 
Air 


Lens 1 


165.470 


9.000 


130.467 





BK7 




98.600 


20.778 


121.809 





Air 


Lens 2 


239.650 


21.356 


121.570 





BK7 




-408.100 


51.233 


119.782 





Air 


Lens 3 


124.500 


7.000 


89.810 





BPM4 




77.790 


15.419 


83.550 





Air 


T ,ATI Q A. 




94 ^"^^ 


81 705 


n 


FPT,51 




-108.870 


4.522 


77.973 





Air 


Lens 5 


-149.580 


6.000 


69.911 





UBK7 




198.530 


17.030 


64.277 





Air 


Filter 


oo 


5.000 


57.061 





BK7 




oo 


10.775 


55.433 





Air 


Camera window 


oo 


4.000 


50.037 





Silica 




oo 


15.329 


48.681 





Vacuum 


CCD detector 


oo 




41.070 
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Fig. 1. — Layout of the four TAOS telescopes at Lu-Lin Observatory. 
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Fig. 2. — Optical layout of the TAOS telescopes. 
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Each telescope is controlled by an Oregon Micro Systems PC68 card. The HA and Dec 
axes are driven by friction wheels which are powered by Bearing Engineers Inc. Advanced 
Vector Servo systems. The axis positions are monitored by Renishaw RGH24 rotary encoders 
with 0.27 arc sec/step resolution. The secondary focus position is set by a screw drive 
powered by an Oriental Vexta PK245-03B motor controlled by an Intelligent Motion Systems 
IB462 bipolar stepping motor driver. 

All of the telescopes were delivered with strong spherical aberrations over the whole 
field as well as significant astigmatism in the corners of the images. The primary mirror 
on the fourth telescope was found to have very severe surface defects which distorted the 
point spread function beyond any reasonable shape and size. Each of the telescopes was 
also found to have severe mechanical and defects. We found the mirrors would move by 
large amounts (the primary mirrors would move 4-5 mm as the telescope would slew from 
horizon to horizon, while the specified tolerance was < 50 fim). We also found the tracking 
accuracy limited by backlash in the drive system reduction gears and bad bearings in the 
friction wheel assemblies. The secondary focus systems were found to have poor resolution 
in positioning and significant backlash. 

We have made significant improvements to the telescopes such that they are now of 
sufficient quality for the requirements of the survey. First, in order to keep the mirrors accu- 
rately positioned, the entire mirror support structures for both the primary and secondary 
mirrors were redesigned and refabricated by engineers at Lawrence Livermore National Lab- 
oratory. Second, the friction wheel assemblies were redesigned and refabricated by engineers 
at the Institute for Astronomy and Astrophysics at Academia Sinica in Taiwan, and the 
reduction gears were replaced with Bayside model PX-34-050-LB (50:1) devices. Third, the 
secondary focus drive systems were completely redesigned by engineers at Academia Sinica. 
A Renishaw RGH22 linear encoder with 0.5 fim resolution was added to each telescope to 
enable the secondary focus positions to be accurately set to within 1.25 fim (the limiting res- 
olution from the stepper motor and threading of the drive shaft). The encoders are equipped 
with limit switches to keep the secondary mirror from bottoming out against the support 
structure in one direction and completely unscrewing from the support structure in the other 
direction. The encoder system is also equipped with a magnetic reference mark to reproduce 
the home position within one motor step, which is needed for the automatic focus control 
(see §5.4.71) . Finally, the primary mirror of the fourth telescope was repolished and resur- 
faced by the Space Optics Group at the Korea Research Institute of Standards and Science 
in Daejeon, South Korea, and it is now of quality comparable to the three other primaries. 

Due to the poor optical quality of the optical components, the actual enclosed energy 
within one pixel averages about 30% rather than the design goal of 60%. The best PSFs 
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now range from 1.6 pixels (4.6 arc sec) in the center of the image to 2.5 pixels (7.3 arc 
sec) in the corners. This lowered the limiting magnitude of the survey. We have partially 
compensated for the reduced optical quality by an upgrade to the CCD camera from a low 
quantum efficiency, front illuminated device to a camera with a back illuminated CCD (as 
will be described in the following section). The limiting magnitude of the survey is now 
roughly R = 13.5, rather than our target of R = 15. We have made a careful search for 
the densest fields along the ecplitic (see fj6]), and we are typically able to monitor fields with 
300 - 500 stars, which is within a factor of 2 - 3 of our design goal. We have designed a 
Hartmann-Shack device to routinely check on the optical performance of telescopes. 



Each telescope is equipped with a thermoelectrically cooled Spectral Instruments Series 
800 CCD camera employing a thinned, backside-illuminated e2v 2048 x 2052 CCD42-40 



chipo The CCDs have 13.5 /im pixels, corresponding to a plate scale of 2.9 arc sec per 
pixel, and the total image area is 27.6 mm x 27.6 mm. There are two read channels, which 
work at a combined speed of 1.8 x 10^ pixels per second. The gain of each amplifier is set 
to about 2.0 ADU/e~, and the read noise is about 10 e~. The cameras are operated at a 
temperature of -20° C, and the dark current is about 0.05 e~/sec/pixel. Each CCD chip has 
the e2v mid-band coating, providing peak quantum efficiency of more than 93% and mean 
quantum efficiency over the spectral range we will employ of over 90%. Each camera has 
a custom-made filter (produced by Custom Scientific, Inc.) which provides more than 80% 
transmission between 500 nm < A < 700 nm. The quantum efficiency of the CCD and filter 
transmission curves are shown in FigureHl 

The camera was delivered with a chiller for the thermoelectric cooler, but we found this 
difficult to operate remotely. Since we do not need chilled water for our application the 
chillers were replaced with ITW Cooling Systems 2500SS welding coolers, which are simply 
water circulators with fan-cooled radiators. Each cooler is equipped with a GEMS 2659 fiow 
switch relay, which will signal if the coolant fiow stops for some reason. 



■^These chips are specified as 2048x2048, but there are four extra rows at the end of the CCD to act as 
a buffer for current leakage from the substrate. These rows must be read out in zipper mode as they are 
shifted with every other row. However, the current leakage is negligible at our 5 Hz image cadence, so we 
can use these extra rows as actual imaging area. 



3.2. 



The Cameras 
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3.3. Weather Monitoring 

Accurate weather monitoring is an important component in any robotic telescope sys- 
tem, especially at a site prone to rain and high humidity. Most of the weather information is 
monitored by two Vaisala WXT510 weather stations. Each weather station is controlled and 
read through a serial port interface by a small Linux workstation. The two systems are used 
for redundancy in case one of them fails. These weather stations provide information such 
as temperature, humidity, dew point, barometric pressure, and wind speed and direction. 

To provide accurate measurement of the dew point under high humidity conditions, we 
have installed a pair of Vaisala HMT337 humidity and temperature sensors. These devices 
are designed for high humidity applications, and specify an accuracy of ±1.7% in relative 
humidity (or more importantly, ±0.3° C accuracy in the dew point). These devices are run 
by the same computers that control the weather stations. 

Each telescope enclosure is also equipped with a Vaisala DRDllA Rain Detector to check 
for any precipitation. The sensors output an analog voltage proportional to the amount of 
moisture on the sensor surface, and these voltages are monitored by a computer in each en- 
closure by a serial analog to digital converter. The devices are extremely sensitive, and will 
trigger the closure of our enclosure lids immediately after the detection of a single droplet. 
Temperatures of the various parts of each optical assembly (primary mirror, secondary mir- 
ror, support struts, etc.) are monitored by Picotech Technology Ltd. thermocouple sensors. 

The weather parameters of most concern are precipitation, humidity, wind, and tem- 
perature. We would like the lids closed if there is any precipitation or if the humidity is high 
enough for condensation to form on the telescope components. Wind is primarily a concern 
because we do not want the lids open if the wind speed is high enough that the telescopes 
might get damaged, or if we are not able to close the lid at all (note that the lid effectively 
becomes a large sail during the opening and closing operations due to the clamshell design, 
described in §3.4p . Wind is also a concern because of the effect on the tracking accuracy 
of the telescopes. As discussed in §11 oscillations in the tracking accuracy due to wind 
will decrease the signal-to-noise ratio achievable by the system. Temperature monitoring is 
necessary because many of the telescope components (e.g. motor controllers) are not rated 
below certain temperatures. Temperatures of the optical and telescope assembly are com- 
pared to the dew point temperature so we can close our enclosure lids if the telescopes get 
cold enough for condensation to form, as well as to provide real time adjustments to our 
focus positions as the temperature changes. (This will be discussed in detail in §5.4.71 ) 

We have therefore defined a set of alarm thresholds for these various weather parameters. 
If one of these parameters crosses a threshold, a 'WEATH_BAD" alarm state is entered. 
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If this alarm state is entered, the hds of all of the enclosures will immediately be closed 
automatically. In addition, we have defined a second set of thresholds which need to be 
crossed for the WEATH_BAD alarm state to be turned off. If a WEATH_BAD alarm state 
is entered, the weather parameters must be past the threshold values for 3 minutes before 
the alarm is turned off. This prevents the alarm from rapidly turning on and off when one 
of the parameters is hovering near the threshold. These thresholds are shown in Table [2l 



3.4. The Enclosures 

The TAOS enclosures (see Figure|5]) were designed to withstand the hurricane force 
winds which accompany the four or five typhoons which typically move through the region 
every summer. The enclosure walls consist of two 1 mm thick steel plates separated by 5 cm of 
insulating foam and are supported by 10 cm x 10 cm steel beams. The working area is about 
3 m X 3 m. The lids are of a counterweighted clamshell type design. They are constructed of 
polyurethane-filled fiberglass and were fabricated by the Aeronautical Research Laboratory 
of the Chungshan Institute of Science and Technology in Taiwan. Each lid is actuated by an 
ADLEE BM-370 brushless DC motor and a 1:400 worm reducer made by Shin Wei Ent Co. 
A custom control circuit is used to provide both manual operation and automated robotic 
operation through software. The time to open or close the lid is about one minute. In order 
to keep the lid operational during a power failure, the motor is powered by a SMR MCS1800 
DC UPS made by Delta Electronics Inc. 

Each enclosure is equipped with air-conditioning to keep the temperature of the telescope 
close to the anticipated night time temperature. This improves image quality by minimizing 
the temperature difference between the telescope components and the ambient temperature, 
and also minimizes the time used to find the best focus and focus change during the night. 

Finally, each enclosure is equipped with a custom isolation transformer to protect the 
hardware from power surges induced by lightning. All four TAOS enclosures are connected 
to the control building by 1 Gbps fiber optic cable for lightning protection as well. 



3.5. The Robotic Observatory Control System 

In each enclosure, every subsystem other than the telescope and camera is controlled 
by a Sealevel Systems Model 8010 32 bit I/O card which is connected to a custom built 



interface box. These boxes are based on a design used by the ROTSE Survey (lAkerlof et al. 



20031 ). but the design has been extensively modified. The boxes control the power to the 
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cameras, the camera coolers, and the telescopes. The boxes also monitor the limit switches 
on the telescopes, and will automatically cut the power to a telescope if one of the switches 
is tripped. This is done in hardware with a series of relay switches to ensure reliability. 
The chiller flow switches are also monitored by these systems, and the software will turn 
off the thermoelectric cooler on the CCD and cut the power to the chiller should the flow 
unexpectedly stop. 

The lids are also controlled by these interface boxes. The lids are opened and closed by 
asserting 12 V signal lines to the lid controller boxes. The limit switches to the lids are also 
monitored, and once again a set of relay switches will automatically cut the power to the 
open/close line when the open/close limit switch is tripped. 

Finally, each box is equipped with an ICS Advent WDT5 watchdog module. This 
module must be pinged by the I/O card once every second, or it will close a set of relay 
switches which will automatically cut the power to the telescope and close the hd. If the 
computer crashes or one of the control daemons should happen to die, the system will then 
be safe from any inclement weather which should arise before the system can be restarted. 

4. High-Speed Photometry with Zipper Mode Readout 

The CCDs used for the TAOS project are described in §3.2[ The devices are 2048x2052 
e2v CCD42-40 chips. Each chip has two amplifiers, which allow for a combined readout rate 
of 1.8 MHz. It takes approximately 2.3 seconds to read out the entire frame, which makes 
our desired 5 Hz image cadence impossible without more than 90% dead time. To overcome 
this limitation, we have designed a novel CCD readout method, which we call zipper mode 
(see Figure[6]). 

As with conventional operation, when reading out in zipper mode, the telescope tracks 
the sky and stars are imaged onto the focal plane. However, instead of simply opening and 
closing the shutter for each exposure, the shutter is left open, and after 105 ms, a "row 
block" of 76 rows is read out. (The read out operation takes approximately 95 ms, so the 
readout time plus the exposure time is equal to our required sampling time of 200 ms.) This 
is followed by another "hold" of 105 ms, then another row block is read out. This cycle 
may be repeated for extended periods of time, possibly up to several hours. Note that after 
27 row blocks are read out (27 row blocks x 76 rows = 2052 rows, the total number of rows 
in the CCD), photoelectrons from every star on the focal plane will be read out in each 
subsequent row block. 

The ability provided by zipper mode to read at 5 Hz comes with a cost in signal-to- 
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noise. This is due to two properties inherent to zipper mode, and any attempt to alleviate 
one of the problems exacerbates the other. First, photons from sky background are collected 
continuously during each hold and shift, while photons from a star are collected only during 
a single hold. Note that with row blocks of 76 rows, the pixels in a single row block contain 
photoelectrons collected from a total of 2052/76 = 27 holds. With a hold time of 105 ms 
and a sampling cadence of 5 Hz, the average recorded number of sky photoelectrons per 
pixel is thus a factor of 27 x 200/105 = 51.4 times brighter than it would be after a simple, 
105 ms exposure. Second, it takes just over 1 ms to read a single row. Photons from stars are 
still collected during this time, meaning that stars will leave streaks as a row block is read 
out, as can be seen in Figure[71 Faint stars are particularly vulnerable to this effect of flux 
contamination by streaks of bright stars. One could conceivably reduce the sky background 
using row blocks with a larger number of rows, but this would reduce the effective exposure 
time relative to the readout time and reduce the measured stellar brightness relative to the 
brightnesses of the streaks. 

We would ideally like the row block size to be an integral factor of 2052, the total 
number of rows on the CCD. (Stars near the edges of an image are missing a fraction of their 
flux, increasing uncertainty on the photometry not only of the edge stars but also of their 
neighbors. This effect is minimized if the row block size is an integral factor of the number 
of rows in the CCD.) We found the signal-to- noise ratio was highest when the read time 
was close to the exposure time, and the integral factor of 2052 meeting this requirement is 
76 rows. 

Because of the increased sky background, zipper mode works well only for relatively 
bright stars. On a dark night, adequate signal-to-noise (~7) can be achieved with a star of 
V ~13.5, which consequently is the faintest magnitude we use in our occultation search. 

This readout mode also requires accurate pointing and tracking throughout the duration 
of a zipper mode run, ideally to within a small fraction (~ 0.1) of a pixel. As shown in 
FigurelHl if there are significant oscillations in pointing error, then neighboring stars in a 
row block can move in different directions from their nominal positions over time. This can 
cause a significant reduction in signal-to-noise for closely neighboring stars in a row block 
due to varying blending fractions over time. 

We did find both significant drifts and oscillations in the tracking of the telescopes, as 
shown in Figure[9l Reducing the magnitudes of these errors required significant hardware 
improvements (see §3.11) and software modifications to the telescopes (see §5.51) . These 
modifications were completed in September of 2005, so a large fr a ction o f our first dataset 



suffered from this problem. This will be discussed in IZhang et al.l (l2008bl ) . 
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5. The Software 

The control software for the TA OS system is based on the control software used for the 



ROTSE Survey (lAkerlof et al.ll2003l ). however it has been extensively modified for TAOS. 
Each component of the system is controlled by a separate daemon. A schematic of the 
daemons and their communication paths can be seen in FigurefTUl and a description of each 
daemon is given in the following subsections. 

Each enclosure has a dual-CPU Pentium 4 computer to control all of the hardware 
described in the previous section, and a small PC104 computer to control the telescope 
via the OMS PC68 Motion Controller Card. The main control computer is located in the 
TAOS control building, which also houses two small computers to monitor the two sets of 
weather stations (see §3.3p . All of the computers are running the Debian Sarg^ distribution 
of Linux. Each of the main enclosure computers are running Linux kernel version 2.4.18 
with the Kansas University Real-Time patch (KURT@) to provide microsecond timing reso- 
lution, which is important for synchronized zipper mode imaging at 5 Hz. The main control 
computer, weather station computers and telescope control computers are running kernel 
version 2.6.8. 

In order to provide accurate timing to the control computers (necessary for accurate 
telescope pointing and synchronized high-speed imaging), the clocks of all of the computers 
in the control system are synchronized via nt;^ to the control computer. The ntp server on 
control computer synchronizes to two ntp servers in Taiwarjll. 



5.1. Command and Status Message Passing 

All command and status information is passed from computer to computer via standard 
inet sockets using a custom message passing library. Command information is sent from the 
control computer directly to each enclosure computer. 

Each computer has a System V shared memory area which contains all status informa- 
tion from each enclosure, the control computer, and the weather computers. We call the 
shared memory area the bulletin board, and the status information is kept up to date by 



|http : //www ■ debian . org 



http : //www . ittc . ku . edu/kurt/ 



^Network Time Protocol. See http : // www . ntp . org/[ 



7 



tick. stdtime. gov. tw, stdtime.sinica.edu.tw 
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the bulletin board daemon, bbd. There is an instance of bbd running on each computer in 
the local network, and once every second each instance of bbd updates every other instance 
with its local status information. Each daemon on the system thus has available up-to-date 
information on every other daemon running on the system. 

5.2. The Scheduler Daemon braind 

The scheduler has four internal states, described in Table [31 While the Sun is above the 
horizon, the scheduler is in the SHUTDOWN state. In this state the lid is closed and the 
telescope and camera are powered off. When the Sun goes below the horizon and the weather 
is good (see §3.3p . the scheduler will enter the GRBMODE state. In this state, the lid will 
be opened and the telescope and camera will be powered on. No zipper mode occultation 
survey data will be collected in this state, however the system will respond to any GRB alerts 
that come in. (Zipper mode data are not so useful in twilight conditions due to the increase 
in sky background, however GRB follow-up data could be useful under this circumstance.) 
When the Sun goes below an elevation of —18° (i.e. the sky is sufficiently dark), braind 
enters STARTUP mode, and it will subsequently commence zipper mode observations for 
the occultation survey. If any GRB alerts come in during a zipper mode run, braind will 
interrupt the observations and schedule follow-up observations for the GRB event. At the 
end of the evening, when the Sun rises above —18°, braind will go back into GRBMODE and 
will then only respond to GRB alerts. At sunrise, braind will go back into SHUTDOWN 
mode for the day, and will close the lids and power down the telescopes and cameras. 

If braind is in the GRBMODE or STARTUP states and at any time the weather should 
turn bad, the system will go into STANDBY mode and close the lids. If the weather becomes 
good again, braind will go into either the GRBMODE or STARTUP state, depending on the 
elevation of the Sun. If the weather is bad when the Sun sets, the system will go straight 
into the STANDBY state from the SHUTDOWN state. Conversely, if the system is in the 
STANDBY state at the end of the night when the sun rises, it will go into the SHUTDOWN 
state. 

Every evening after sunrise and sunset, braind will schedule a set of dark images to 
be taken on each camera. At sunset, this happens 300 seconds after braind enters the 
GRBMODE state to give the camera ample time to cool down. At sunrise, the dark images 
are taken after the lids are closed but before the system enters the SHUTDOWN state. Each 
series of dark images comprises a set of 20 stare mode dark images with 1 second exposures, 
and a series of 10 zipper mode dark images with the standard zipper mode parameters (5 Hz 
sampling, blocks of 76 rows, and 32 row blocks per FITS file). 
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Table 2. Weather alarm parameters 



PciT'CiTYiO'J'OT' A 1 QyTTi (in T^lrroanr^lH 
JT dl dlliCtd jT-idllli vyil ± 111 tJbllUiLl 


S\LCXL 111 ii A. Ill CollUiLl 


DRDllA Voltage < 2.5 V 


> 2.8 V 


Relative Humidity*^ > 95% 


< 94% 


Dew point - temperature differential^ < 1.0° C 


> 1.5° C 


Average wind speed > 20 km/hr 


< 18 km/hr 


Maximum wind speed'^ > 30 km/hr 


< 28 km/hr 


Minimum temperature <-5.0° C 


> -3.0° C 


^Value reported by HMT337 dew point sensor. 




''Dew point value reported by HMT337 dew point sensor, 


temperature differentials 


calculated from HMT337 temperature reading and Picotech thermocouple measure- 


ments of telescope component temperatures. 




^Maximum wind speed measured over a 10 minute interval. 




Table 3. Scheduler states 




State Sun Elevation Weather Status Lid Position 


Telescope and Camera Power 


SHUTDOWN > 0° good or bad closed 


off 


GRBMODE < 0° good open 


on 


STARTUP < -18° good open 


on 


STANDBY <0° bad closed 


on 
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5.3. The Weather Daemon bbweathd 

The weather stations and dew point sensors described in §3.31 are read by two instances 
of bbweathd, one instance for each set of one weather station and one dew point sensor. This 
daemon (as the name imphes) simply reads both devices and pubhshes the information on 
the bulletin board for other software components to monitor. In addition, bbweathd will 
check the measured weather parameters against a set of limits, and set an alarm if adverse 
weather is detected. This alarm state is also written to the bulletin board for other daemons 
to monitor. 



5.4. Control Daemons 

In addition to the instance of bbd described above, each enclosure computer runs a 
set of eight control daemons, each dedicated to a particular hardware component (or set 
of components) or high level control task. These daemons are described in the following 
subsections. 



5.4- 1- The Master Control Daemon taosd 

Every daemon running on the enclosure computer is a child of taosd. This daemon is 
responsible for forwarding commands from braind to the appropriate child daemon, and for 
monitoring the status of each of the children and passing the status information on to the 
local instance of bbd. Each child daemon has a dedicated System V shared memory area 
which is used to pass command and status information to and from taosd. This daemon will 
automatically send commands to certain daemons based on the status information of other 
daemons (e.g. it will close the lid if adverse weather is detected). 

5.4-2. The Watchdog Daemon spotd 

This simple daemon is responsible for pinging the watchdog module in the observatory 
control interface box described in §3.5[ The daemon is designed to run independently of the 
rest of the daemons on the system to enable manual control (e.g. so the lid can be opened 
manually without the rest of the software running), however when taosd is running, a simple 
"I'M ALIVE" command is expected from taosd at a rate of 1 Hz. If this command times out, 
spotd will die, the telescope will be stopped and the lid will automatically be closed. Thus 
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if taosd should crash or otherwise be hung up for some reason, or if the enclosure computer 
itself should crash, the system will automatically go into a safe mode. 

5.4.3. The Local Weather Daemon weathd 

While general weather information is monitored by bbweathd as described above, a local 
weather daemon weathd runs on each enclosure computer as well. This daemon has three 
tasks. First, weathd monitors the Vaisala precipitation sensors described in §3.31 Second, it 
monitors the local copy of the bbweathd information on the bulletin board. If precipitation 
is detected or bbweathd reports adverse weather, a flag is set in the weathd shared mem- 
ory segment which will instruct taosd to abort any current observations and close the lids. 
Finally, weathd monitors the telescope strut and mirror temperatures read by the thermo- 
couple sensors attached to the telescope (see §3.ip . and compares these temperatures to the 
dew point reported by bbweathd to see if the lid needs to be closed to prevent condensation 
on the telescope. 

5.4 ■4- The Telescope Command and Status Daemon torusd 

The telescope is directly controlled by a small PC104 computer via an OMS PC68 motion 
controller card, and communication between the main enclosure computer and the PC104 
computer is handled by torusd. This daemon connects to the telescope control daemon teld 
(see §5.51) over a set of four inet sockets, one for pointing and tracking control, one for the 
focus control, one for status information and one for miscellaneous control commands. 

5.4-5. The Guide Daemon guided 

While we have greatly improved the pointing and tracking capabilities of the system 
through both software and hardware modifications, we found the pointing was still inade- 
quate for the survey. Accurate pointing (within 0.3 arc sec) must be maintained throughout 
an entire zipper mode run, which can last as long as 2 hours. As discussed in 21 poor track- 
ing can significantly decrease the signal-to-noise performance of the system, so it is desirable 
to keep the tracking as accurate as possible. While many surveys have a separate guiding 
imager, TAOS has the advantage that images are collected at a high rate, and we use this 
feature to enable our guiding system. 

The guide daemon runs in two modes. During the initial pointing phase, we point the 
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telescope to a field center and take a stare mode image. This image is analyzed to find 
the actual RA and Dec of the image center, and a pointing offset is calculated. A pointing 
adjustment command is sent to torusd, and the process is repeated until the center of the 
image is within 0.1 pixels of the actual field center. In this way, at the beginning of a zipper 
run, all four telescopes are pointing to within 0.3 arc sec of the target field center. 

The image analysis routine makes use of the standard software packages SExtractor 
(IBertin fc Arnoutsi Il996r) and w cstools (IMinkl |2006| ) used in conjunction with the USNO- 
Bl.O catalog (IMonet et al.ll2003l ). First, after an image is taken, the wcstools utility imwcs is 
run to find the World Coordinate System (WCS) keywords to transform RA and Dec to pixel 
coordinates. Second, SExtractor is run to find objects in the image. Third, the SExtractor 
output is used with the wcsto ols utility immatch to mat ch the objects found by SExtractor. 
Finally, a set of polynomials (iMalumuth fc Bowers! 119971 ) are fit to the data to convert from 
X and y, the pixel coordinates of each matched star calculated from the WCS coefficients, to 
Xp and yp, the actual centroid locations reported by SExtractor. The polynomials are of the 
form 



Xr. 



Co + Cix + C2y + C^x^ + dxy + C^y^ 
+C6X^ + Crx^y + Cgxy^ + Cgy^ 
Do + Dix + D2y + D^x^ + D^xy + D^y'^ 
+Dex^ + Djx^y + D^^xy'^ + Dg?/^ 



where and Di are the fit coefficients. After the fits are complete, these distortion coeffi- 
cients and WCS keywords are used to calculate the RA and Dec of the center of the image, 
which is then used to calculate the pointing offset. The pointing offset is written to the 
guided shared memory status segment, which is then read by taosd and passed to torusd as 
a pointing adjustment command. 

The second guided mode is the tracking mode. At the end of the initial pointing phase, 
20 guide stars are selected from the stare mode image. The stars are selected on criteria 
including moderate brightness (i.e. bright enough to have good signal-to-noise characteristics 
but not so bright as to leave significant streaks in the zipper images), isolated (no close 
neighbors in zipper mode images), and they should be spread out through the image so the 
entire focal plane is well sampled. Note that a set of 32 consecutive row blocks is typically 
written to each FITS file, so there are effectively 640 guide stars in each image. Every 
30 seconds, the last FITS image that was written to disk is analyzed. SExtractor is run on 
the image, and the centroids of the guide stars are found and compared with the original 
values. A pointing correction is calculated and written to the guided shared memory status 
segment, which in turn is read by taosd and passed as a pointing adjustment command to 
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torusd. We find that the inherent tracking of the system is sufficient that the pointing drifts 
by much less than 0.1 pixel in the 30 second guiding update period, and guided is thus able 
to keep the pointing correct throughout an entire zipper mode run of up to two hours in 
duration. 

The guide daemon also monitors image rotation throughout the course of a zipper run. 
Taiwan suffers from frequent earthquakes which can disrupt the polar alignment of the 
telescopes. We can thus learn in real time if the telescopes need realignment. Furthermore, 
guided reports the visibility conditions throughout the duration of a zipper mode run by 
monitoring the star count. Each time a zipper mode FITS file is analyzed for tracking 
adjustments, SExtractor reports the number of stars to which it can fit a standard PSF. The 
star count is compared to the number of stars found at the beginning of the zipper mode run 
to indicate how the observing conditions vary throughout the run. In our fully automated 
observing mode, taosd will interrupt observations if the star count drops significantly. This 
usually occurs due to incoming clouds, or if the seeing degenerates due to humidity or 
focusing problems. This real-time image quality analysis helps optimize our data collection 
as we can interrupt a zipper mode run and refocus the telescope if necessary. 

5.4-6. The Lid Control Daemon clamd 

The lid is opened and closed by a relatively simple daemon clamd. This daemon will 
assert an open or close signal to the lid controller on command, and monitor the open and 
close limit switches on the lid. The open and close signals are lowered when the appropriate 
limit switch is hit. The daemon will also set an alarm if the open or close operation times 
out, and this alarm will alert someone to fix the problem manually. 

5.4-7. The Focus Control Daemon focusd 

The focus control daemon focusd has two tasks. The first is to analyze images taken 
during a focus sequence to determine the ideal focus position. A focus sequence consists of 
a series of one second exposures taken with secondary mirror at different focus positions. A 
total of 17 images are acquired while moving the focus position by 15 fim. steps. A Ikxlk 
sub-image from the center of the focal plane is analyzed using SExtractor to find the average 
Pull Width Half Maximum (FWHM) of each star in the image, and the focus position 
and medium FWHM are logged internally. After all of the images have been collected and 
analyzed, the FWHM data are fit to two straight fines, one line on each side of the point with 
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the minimum FWHM. These fits are then used to calculate the best focus position, which 
is then entered into a database along with the telescope component temperatures (reported 
by weathd) and zenith angle of the telescope (reported by torusd). 

The second function of the focus daemon is to monitor the telescope component tem- 
peratures and zenith angle throughout a zipper mode run and to look up the ideal focus 
position in the database. If the ideal focus position changes, the new position is written to 
the focusd shared memory status segment, where it is read by taosd and passed on to torusd 

focus adjustment command. 



5.5. The Telescope Control Daemon teld 

The system originally used the telescoped program from the OCAA^ software package 
to control the telescope. However, we found three problems with this code that led us to 
the decision to write a completely new software package. First, the pointing model used by 
OCAAS did not accurately represent the pointing actually realized with the system. The 
residuals between the model and the true pointing were often as high as 10 arc minutes. 
Second, we found it difficult to generate the pointing model under OCAAS, where the point- 
ing is adjusted by hand until a catalog star is centered on the image, and this is repeated 
for several stars. This was very time-consuming as it needs to be done several times per 
year (after earthquakes), and it needs to be done for all four telescopes. It is also difficult 
to do remotely. Third, we found it difficult to implement a robust communication interface 
between torusd and the OCAAS software without a significant rewrite of the software. 

The replacement daemon, teld is a multi-threaded program with separate threads dedi- 
cated to telescope pointing, focus control, status reporting, and master control (i.e. a parent 
thread to monitor the children). In addition, when tracking a field a new thread is spawned 
to control the telescope tracking directly. This thread runs with an elevated scheduling prior- 
ity to keep the timing of the control loop accurate. Implementing a multi-threaded solution 
to the telescope control also required rewriting the driver for the PC68 motion controller 
card to make it thread safe. This was done by simply adding a semaphore to the driver 
which is raised whenever the card is accessed. 

The daemon makes use of the SLALIB Positional Astronomy Librar}]^ to convert to and 
from J2000 RA and Dec to observed HA and Dec (/lobs and 5obs)- A custom pointing model 



® Observatory Control and Astronomical Analysis System, see 
* ^http : //star-www . rl . ac . uk/ 



http : //www . clearskyinstitute . com/ 
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derived from that used by the TPOINT software packag is used to convert hohs and 5obs 
to actual encoder angles {6h for the hour angle encoder and 6s for the declination encoder). 
The pointing model used is: 

hi = Khs - TF cos 6*1^1 sin hohs/ cos 5obs 
5i = f^obs - TF(cos 6ii,t cos hohs sin 6ohs 

-sin 6'iatcos5obs) 

— KZ cos 5obs — MZ sin 6ohs 

-FO cos hohs 
= -hi + CH sec 6i + NP tan 6i + 



9, 



(ME sin hi - MA cos hi) tan 5i + IH 
= 5i - ME cos hi - MA sin /ii - ID 



where hi and 6i are temporary variables, ^lat is the geodetic latitude of the site and the ten 
pointing model coefficients are described in Table |H 

After implementing teld, we found the typical residuals between the calculated and true 
positions to be less than 10 arc seconds. Used in conjunction with guided, we found that we 
can track to well within 0.3 arc sec for durations up to at least two hours. 

We also implemented an automatic pointing model generation feature to the system. A 
set of 40 images is taken at a grid of locations on the sky, and the images are analyzed by 
guided to find the actual RA and Dec of the centers of the images. The subseque nt arrays of 



RA, D ec, 9h, Os, and epoch are fit using the amoeba fitting routine adapted from lPress et al. 



(I1994J ). This works much faster than the OCAAS method, especially when it must be done 
to all four telescopes several times a year. This also allows a new pointing model to be 
created remotely and automatically. 



5.6. Collection of Zipper Mode Data 

When the system is in STARTUP mode, braind will schedule a series of zipper mode 
runs. At the start of a run, braind will choose a standard TAOS field (see based on 
criteria such as distance from the ecliptic plane, distance from zenith, and number of stars 
in the field. The scheduler will then instruct torusd to point all of the telescopes at the field 
center. Each telescope will then undergo a focus sequence (see §5.4.71) and set the secondary 



http : //www . tpsof t . demon . co . uk/ 
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focus to the best position. Once this is completed, each telescope will undergo a pointing 
correction, as described in §5. 4. 51 Once all of the telescopes are pointing correctly, a set of 
three 1 second stare mode images are acquired for diagnostic purposes. Next, the zipper 
mode sequence is started. Each camera system is sent a start time by braind, and when the 
start time is reached, the shutters are opened and zipper mode readout is commenced. This 
5 Hz readout process continues for a total of 90 minutes (as long as the weather remains 
good and no GRB alerts are received). Throughout the duration of the run, guided will keep 
the telescope pointed correctly. At the end of the run, a second set of three diagnostic stare 
mode images are acquired, also with 1 second exposures. Finally, a sequence of 1 second 
stare mode images is acquired with pointing offsets of 0.8° on a 3x3 grid about the center 
of the field. This provides an effective mosaic image of ~lin° centered on the field, which 
is also used for diagnostic purposes. 

Throughout a zipper mode run, braind will monitor the results of each command and 
the status of each daemon. If an error occurs, braind will attempt to recover by restarting 
any software that needs it and, if necessary, reinitializing the zipper mode run. If guided 
reports bad observing conditions, the zipper run will be aborted. A new field may be selected 
at this time based on the scheduling criteria discussed above, and if so braind will instruct 
the telescopes to point to the new field. In any case, the system will then undergo a new 
focus sequence to refocus the telescopes. If the weather is cloudy, this step will fail due 
to the lack of stars in the images. In this case, the system will wait for 300 seconds and 
reattempt the focus sequence. This step will be repeated until the weather clears. If the 
cloudy weather lasts long enough, a new field may be selected by braind in the meantime. 

Every morning, all of the FITS images from the previous night (typically about 40 GB 
per telescope on a clear night) are archived to a local RAID server in the main control 
building. Once every month or so, the data are copied to a set of hard drives which are 
physically brought to Academia Sinica. The data are copied to a local data server, and then 
transfered to the Taiwan National Center for High- Performance ComputinJ^ for backup to 
tape. Finally, the data are copied over the internet to a second RAID server located at 
National Central University. We then have three copies of the data in separate locations for 
secure, redundant data storage. The data are subsequently analyzed on using two computer 
clusters (at Academia Sinica an d National Central U niversity) with a custom photometry 



package that will be described in IZhang et al.l (l2008bl ). 



http : //www . nchc . org . tw/ 
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5.7. The GRB Alert Response System 

The TAOS GRB response system consists of two components, the forwarder and the 
GRB daemon grbd. The forwarder simply listens to a socket connected to the GCN, and 
passes the information on to grbd. The purpose of the forwarder is to run on a machine 
with a more reliable internet connection than that available at Lu-Lin Observatory. If the 
internet connection to Lu-Lin is down, the forwarder still receives the GCN packet and keeps 
trying to forward the packet to grbd until it succeeds. We currently run the forwarder on a 
computer at the Institute for Astronomy and Astrophysics at Academia Sinica in Taipei. 

When grbd receives a GCN alert packet from the forwarder, it parses the message and 
determines whether the reported event can be followed up by the TAOS system. If the GRB 
event is observable, grbd calculates an observation schedule and passes the request on to 
braind, which will interrupt any zipper mode operations and schedule the GRB follow-up 
observations. 

TAOS has the advantage of multiple telescopes operating at the same site, which makes 
it possible to use different exposure times for the same event. Our observing scheme is 
optimized to have at least one telescope looking for fine structure in the light curves of bright 
events with short exposures, and two telescopes taking longer exposures for sensitivity to 
fainter afterglows. In the future, we may schedule one telescope to observe in zipper mode 
at 5 Hz to look for very bright and rapidly varying transients. To date, 10 GRB alerts 
were resp onded to by the T AOS system, and two afterglows (R = 16.5 and 16.8) have been 



detected flWang et al.ll2008h 



6. Field Selection 

In order to meet the survey requirements discussed in Hl.l\ a set of fields needed to be 
found which had sufficient numbers of stars. This task was complicated by the fact that the 
TAOS survey would be using zipper mode readout, as the brighter stars will leave significant 
streaks in the images. It was therefore desirable to minimize the numbers of bright stars in 
the fields. A set of 167 fields was found, with varying characteristics such as number of stars, 
ecliptic longitude and latitude. Information on the TAOS fields can be found on the TAOS 
websit^. Note that this set is much larger than is needed, however a large set was generated 
for the first phase of the survey in order to evaluate the fields in terms of signal-to-noise as 
a function of parameters such as sky background, stellar density, etc. After analysis of the 



^ ^http : //taos . asiaa. sinica. edu . tw/taosf ield/ 
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first data set, a small subset of these fields will be selected for further observations. 



The fields were selected by using the USNO Bl.O catalog (iMonet et al.ll2003l ). Initially, 
lists of stars were generated in 300x300 arc second tiles. These tiles covered all of the sky 
with declinations in the range of —22° < 6 < 68°, which corresponds to elevation angles 
higher than 45° at the Lu-Lin site. After generating over 4.6 million tiles, these smaller 
sections of the sky were merged into larger 1.7° x 1.7° lists of stars. The center of each small 
tile (except for those near the declination limits) became the center of the larger field. This 
process ensured overlap, allowing for the ability to maneuver around the bright stars {R < 9) 
in the sky. Over four million of the large fields were generated. 

A subset of these fields was then selected according to the following characteristics: 



• Maximum number of stars in the field with i? < 15. 

• Minimum number of bright stars {R <9). 

• Maximum coverage of the ecliptic. 

• Some control fields away from the ecliptic. 

• No nebulae or other astronomical phenomena. 



These requirements were used to narrow down the list of viable fields. The field list was 
refined using a sorting routine to select fields with varying star densities. The main selection 
requirement was minimizing the streaks in the zipper mode images while not significantly 
compromising on the total number of stars. Candidate fields that did not meet these criteria 
were cut from the list. Fields with nebulae or other phenomena were removed by checks of 
Digital Sky SurvejQ (DSS) images. The DSS images were also used to verify the accuracy 
of the USNO Bl.O catalog, particularly to ensure that there are no bright stars in the field 
that do not appear in the catalog. 

Note that it was impossible to meet the stellar density requirement in all of the fields, 
as most of the fields along the ecliptic are very sparse (some of these fields have fewer than 
100 stars). However, care was taken to find the densest fields available in these sparse regions, 
and it should be noted that the sparsest fields are visible mostly during the rainy summer 
months when observing is limited. Also note that high priority is given to observations of the 
densest fields in order to maximize the amount of useful data collected, and that the sparser 
fields are only observed when necessary due to the season and phase of the moon. The 



^ ^http : //archive . stsci . edu/dss/] 
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typical number of st ars in the fields observed in our first data set was thus 300 - 500 stars 



fIZhang et al.ll2008af ) 



7. System Performance 

To test the overall performance of the system we have observed several predicted asteroid 
occultation^. Figure [11] shows the lightcurves of the occultation of the star TYC 076200961 
{V = 11.83) by the asteroid Idea {V = 14.0, diameter = 97 km) at UTC 12:35 February 
6, 2006. The event was observed by operating the three functioning telescopes at 4 Hz, 
and the event, which lasted approximately 6 seconds, is clearly visible, indicating that the 
system is indeed capable of accurate, high-speed imaging. Figure[T2] shows a series of zipper 
mode images of the occultation of the star HIP 050535 (V^ = 8.46) by the asteroid Klemola 
{V = 15.7, diameter = 31 km) at UTC 12:10 June 5, 2004. This was a short duration event, 
lasting only 1.3 seconds. 

These asteroid occultation events clearly illustrate that the TAOS system is capable 
of detecting short duration occultation events with multiple telescopes. The system has 
been collecting data for over three years with three telescopes, and the fourth telescope 
came online in August 2008. Since February of 2005, we estimate TAOS has made over 
5 X 10^° individual photometric measur ements, and analysis of the first two years of data 



has been completed (IZhang et al.ll2008al ). 
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Fig. 3.— The TAOS B telescope. 
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Fig. 4. — Spectral response of the camera system. The dashed hne shows the quantum 
efficiency of the CCD, the dotted hne shows the ffiter transmission, and the sohd hne shows 
the total spectral response. The quantum efficiency data for the CCD was obtained from 
the e2v website at http : //www . e2v . com[ 
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Fig. 5. — The enclosure for telescope TAOS A. 
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Fig. 6. — Illustration of zipper mode operation with the focal plane divided into four row 
blocks. Solid circles indicate actual stars, empty circles indicate locations of electrons after 
shifts. Note that after the fourth sequence of expose and shift operations, flux from all four 
stars is read out at each subsequent shift, although the flux from these stars was collected 
at different epochs. 
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Fig. 7. — Example of a zipper mode image comprising a 512 x 76 subsection of a 2048 x 76 
row block. 



Table 4. Pointing model coefficients 



Coefficient Description 



IH 


Hour angle home switch position 


ID 


Declination angle home switch position 


CH 


CoUimation error 


NP 


Non-perpendicularity of h and S axes 


MA 


Azimuthal polar alignment error 


ME 


Elevation polar alignment error 


TF 


Tube flexure 


FO 


Fork flexure 



KZ and MZ First order correction to declination encoder error 
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Fig. 8. — Continuation of zipper mode schematic shown in Figure O illustrating effect of 
traclfing errors throughout a zipper mode run. In this illustration, the telescope drifts to the 
left during exposure 7^8, moving the stars to the right of the image. During exposure #9, 
the telescope moves too far back to the right, moving the stars to the left of the image. The 
arrows on the images indicate the displacements of the stars from their nominal positions. 
Note that the stellar ffux that is read out in a row block is collected at different times, 
depending on where a star is located on the focal plane. Neighboring stars on a row block 
could therefore have offsets from their nominal positions in different directions. 



Fig. 9. — Two row block subsections taken at two different epochs. Note the two pairs of 
stars inside the circles. The two neighboring stars in the right circle are close together in the 
left image and have moved apart in the right (later) image. The two stars in the left circle 
moved closer together in the right image. 
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Fig. 10. — Schematic of the TAOS system control software. Sohd arrows indicate control 
message bus, and dashed arrows indicate status message bus. See text for details. 
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Fig. 11. — Lightcurves from three TAOS telescopes of the occultation of the star 
TYC 076200961 {V = 11.83) by the asteroid Idea {V = 14.0, diameter = 97 km) at UTC 
12:35 February 6, 2006. 
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Fig. 12. — Zipper mode sub- images from telescopes TAOS A and TAOS B of the occultation 
of the star HIP 050535 {V = 8.46) by the asteroid Klemola {V = 15.7, diameter = 31 km) 
at UTC 12:10 June 5, 2004. Each sub-image is an exposure of 0.25 seconds, with time 
increasing to the right. The vertical lines visible in each set of images are caused by a slight 
excess in dark current from the CCD substrate, and indicate the boundaries of the row 
blocks. These images were taken before guided was implemented, so the two telescopes have 
different pointing offsets and the occulted star is in a different location in the row block. 



